Optimization of evaporative cooling towards a large number of Bose-Einstein 

condensed atoms 
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We study the optimization of evaporative cooling in 
trapped bosonic atoms on the basis of quantum kinetic the- 
ory of a Bose gas. The optimized cooling trajectory for *^Rb 
atoms indicates that the acceleration of evaporative cooling 
around the transition point of Bose-Einstein condensation is 
very effective against loss of trapped atoms caused by three- 
body recombination. The number of condensed atoms is 
largely enhanced by the optimization, more than two orders 
of magnitude in our present calculation using relevant experi- 
mental parameters, as compared with the typical value given 
by the conventional evaporative cooling where the frequency 
of radio-frequency magnetic field is swept exponentially. In 
addition to this optimized cooling, it is also shown that highly 
efficient evaporative cooling can be achieved by an initial ex- 
ponential and then a rapid linear sweep of frequency. 

Pacs. 03.75.Fi, 05.20.Dd, 32.80.Pj 



I. INTRODUCTION 

Evaporative cooling is an essential experimental tech- 
nique in the recent demonstrations of Bose-Einstein con- 
densation (BEC) with magnetically trapped atomic va- 
por . This cooling method is based on both the 
selective removal of energetic atoms through evaporation 
and collisional rethermalizations among remaining atoms 
1^,0 . The successful utilization of interatomic elastic col- 
lisions makes evaporative cooling highly efficient. A real 
system, however, involves other atomic collisions, such as 
elastic collisions with background gas, inelastic collisions 
due to dipolar relaxation and three-body recombination 
processes. These undesirable collisions result in a loss of 
the trapped atoms and seriously reduce the effectiveness 
of evaporative cooling. 

For alkali-metal atoms, such as *^Rb and ^■^Na, three- 
body recombination becomes a dominant loss mechanism 
during evaporative cooling and limits the final num- 
ber of condensed atoms, Nq, in the experiments. It is 
known that *^Rb atoms suffer much severer losses than 
^^Na atoms due to their much large rate coefficient of 
three-body recombination process [||J|]. BEC experi- 
ments have clearly reflected this fact: the reported A^o 
value in ^'^Rb is 10'^ at most pOl], while that in ^SNa is 



up to 5 X 10'^ 0. The small BEC of ^^Rb is a great 
disadvantage with respect to its future application as an 
"atom laser" |12 and has to be overcome. 

In this paper, we report the theoretical optimization 
of evaporative cooling, focusing on BEC experiments 
with *^Rb atoms. Several theoretical works have con- 
sidered the optimization of evapor ative cooling to in- 
crease the cooling efficiency |l3yi4| ,|6|. Especially, Ref. 

optimized the evaporative cooling in atomic hydro- 
gen not only theoretically but also experimentally. How- 
ever, these studies are based on kinetic theory of a classi- 
cal gas which fails around the BEC transition point, and 
therefore don't give us any information on the number of 
condensed atoms produced by evaporative cooling. Our 
analysis is based on quantum kinetic theory of a Bose 
gas [|5| and can provide quantitative results on the in- 
crease in the number of condensed atoms by optimizing 
the evaporative cooling. We optimize the temporal fre- 
quency sweep of the radio-frequency (rf) magnetic field 
under a fixed trapping potential by following the vari- 
ational method developed by Sackett et al. for evap- 
orative cooling in a classical gas system [^. We show 
that an accelerated sweep of the rf-field frequency around 
the BEC transition point can largely enhance the final 
number of condensed atoms compared with conventional 
cooling with an exponential sweep of the frequency. It is 
also demonstrated that highly efficient evaporative cool- 
ing can be achieved by an initial exponential and then 
a rapid linear sweep of frequency. Our present results 
provide a useful guideline to experimentalists aiming to 
achieve large condensates in ^^Rb atoms. 



II. KINETIC THEORY OF EVAPORATIVE 
COOLING 

A. Thermodynamics of trapped bosonic atoms 

We begin with a review of our theoretical formulation 
p5[ . During evaporative cooling, a magnetic potential, 
IJjv) , is truncated depending on the frequency of the ap- 
plied rf-magnetic field. The thermalized distribution of 
noncondensed atoms in such a truncated potential is well 
approximated by the truncated Bose-Einstein distribu- 
tion function such that 
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where T is the temperature, tp — 12171 is the kinetic 
energy of atoms with momentum p and mass m, and 
Q{x) is a step function. In this paper, we add tilde to 
the notations of the quantities evaluated through this 
truncated distribution function, ^(r) denotes the local 
fugacity including mean-field potential energy; ^(r) = 
exp{[^ — U{r) — 2vn(r)]/kBT}, where fj, is the chemical 
potential, v = Anah /m is the interaction strength of 
atoms in proportional to the s-wave scattering length a, 
and n(r) is the density profile of atoms. The step function 
eliminates the momentum states whose kinetic energy 
exceeds the effective potential height A{r) = — U{r) ~ 
2vh(r), where et is the truncation energy of a magnetic 
potential. We calculate the density profile of atoms using 
this truncated distribution function in a self-consistent 
way such that 



n(r) =4TTh ^ J f{r,p)p^dp, 
and the internal energy density similarly as 



(2) 



g(r) = ATTh"^ J epf{r,p)p'^dp + vfi^{r) + J7(r)n(r). (3) 

The spatial integrations of these density functions give 
the total number of atoms, N ^ j h{r)dr, and the total 
internal energy, E — J e{r)dr respectively. 

After the BEC transition, the atom density in the con- 
densed region is obtained by the sum of a condensate 
no(r) and saturated noncondensed atoms n„ such that 
n(r) = no(r) -|- n„. The condensate is described by the 
Thomas- Fermi distribution no{r) = Up — U{r)/v with 
the peak density of the condensate at the center of po- 
tential rip, while ri„ is evaluated through the truncated 
Bose-Einstein distribution function with the condition 
|(r) = 1. 



B. Dynamics of evaporative cooling 

The dynamics of trapped atoms during evaporative 
cooling is investigated on the basis of quantum kinetic 
theory of a Bose gas, developed previously by the present 
authors Trapped atoms are removed from the po- 

tential by not only evaporation but also undesirable col- 
lisions, as mentioned in the introduction. The change 
rates (i.e., loss rates) of the total number of atoms and 
the total internal energy are calculated respectively as 
the sum of the contributions of all processes such that, 
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^ = -J "cv(r) dr 

G.|if,(r)[n(r)rdr+(^|^^ e„ 

J eev(r)dr (4) 
-pGsJ if.(r)e~(r)[n(r)]-i + et. 

The rates tt-cv and e^v denote the evaporation rates of 
density functions derived from a general collision inte- 
gral of a Bose gas system ||l5| , [l6|| . The parameters Gi, 
G2, and G3 are the decay rate constants of trapped atoms 
due to the background gas collisions, dipolar relaxation, 
and three-body recombination, respectively. Ks repre- 
sents the correlation function which describes the s-th 
order coherence of trapped atoms, and we assume the 
expressions of Kg for an ideal Bose gas system JT7|p|. 
The terms proportional to the change rate of truncation 
energy, it = det/dt, give the contribution of extra atoms 
"spilled over" when et changes continuously in the forced 
evaporative cooling . The dynamics of the evapo- 

rative cooling including the loss of trapped atoms is de- 
scribed by the coupled differential equations in Eq. (Q). 
In our theoretical framework based on truncated Bose- 
Einstein distribution function, all thermodynamic quan- 
tities are given as the complicated functions of three in- 
dependent variables: temperature T, chemical potential 
/i, and truncation energy e^. We derive here the useful 
thermodynamic relations in terms of the change rates of 
these variables: 
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where T and are the change rate of temperature and 
that of chemical potential respectively. Substituting 
these relations into Eq. , we clearly see that the terms 
proportional to et are canceled out and both T and /i are 
determined by three variables T, /i, and et- 



C. Optimization of cooling trajectory 

The optimum cooling trajectory in the evaporative 
cooling process is obtained approximately by following 
the instantaneous optimization procedure investigated in 
Ref. II^]. We determine the optimized cooling trajectory 
so as to achieve the maximum increase in phase-space 
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density p with the smallest decrease in the number of 
trapped atoms N . In the nonuniform gas system, the 
phase-space density is evaluated at the peak position and 
given by p = fi[0)\^ with the thermal de Broglie wave- 
length A — h/^/2TrmkBT. At each time, we calculate the 
optimum truncation energy by maximizing the deriva- 
tive, —d(\n p) I (i(ln jV) , which describes the efficiency of 
evaporative cooling fl^,^. Wc use a more useful form 
of this derivative rewritten in terms of change rates of 
temperature, chemical potential, and truncation energy 
as, 
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In the optimized evaporative cooling, the truncation pa- 
rameter 77 — et/ksT changes more slowly than the vari- 
ables T, /i, and et- We neglect the change rate of this 
truncation parameter, and consequently et is related to T 
as et/e* = T /T in Eq. (^. The value of optimum trunca- 
tion energy is given as a function of both temperature 
T and chemical potential p.. Substituting this into Eq. 

, we obtain the coupled differential equations describ- 
ing the dynamics of the optimized evaporative cooling. 
The time-evolution calculations are carried out numeri- 
cally by repeating a series of steps consisting of instanta- 
neous optimization and dynamical short-time evolution 
of Eq. (|) with the temporal value of e^. Furthermore 
we assume quick rethermalization, where the system is 
always described by the truncated distribution function 
[p5|,[7[. The obtained function el{t) corresponds to the 
optimized sweeping function of the rf-field frequency. 

In the next section, we will show the numerical re- 
sults calculated by using relevant experimental parame- 
ters and discuss the experimental possibility of making 
large condensates in ^^Rb atoms. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

A. Optimized evaporative cooling 

In our numerical calculations, we employed parame- 
ters corresponding to the experiments performed at the 
University of Tokyo [^. The cloverleaf magnetic trap 
with bias field Bq, radial gradient B , and axial curvature 
B was modeled by the ma gnetic field of loffe-Pr itchard 
type such that _B(r, z) = y^(ar)2~+"(^z^"^r~Bj^, where 
a = - B"Bo/2 and f3 ^ B" /2. We used the typ- 

ical values adopted experimentally, such as Bq = 0.9 G, 
b' = 180 G/cm, and b" = 170 G/cm^. The corre- 
sponding trapping frequencies are ujr = 2t: x 171 Hz in 



the radial direction, and lOz = 27r x 11.8 Hz in the ax- 
ial direction. The coUisional parameters for the trapped 
F = 1, mp = — 1 state were set as follows s-wave scat- 
tering length a ~ 5.8 nm, loss rate constants due to back- 
ground gas collisions Gi = 0.01 s^^, dipolar relaxation 
G2 — 1.5 X 10^^® cm'^s^^, and three-body recombination 
G3 ~ 4.3 X 10^2^ cm®s~^. We chose the temperature 
T = 150 pK, the peak density h{0) = 4.0 x 10" cm'^^ 
and the number of trapped atoms iV = 1.0 x 10^ as the 
initial evaporative cooling conditions. 

The exponential sweep of the rf-ficld frequency is usu- 
ally adopted in BEG experiments ||l^ . We also performed 
calculations assuming that the frequency, u, is swept ex- 
ponentially from 35 to 0.67 MHz in 60 s as in Ref. [p^ . 
The comparison between the different sweeps of rf-field 
frequency, optimized or exponential, will clarify the merit 
of the present optimization procedure. 
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FIG. 1. Sweeps of rf-field frequency in the optimized case 
(solid) and in the exponential case (dashed). The inset shows 
the last part of the optimized sweep using an expanded scale. 
The arrow in the inset indicates the point at which the BEC 
transition occurs. 



Figure 1 shows each sweep of the rf-field frequency. We 
see that the sweeping time (i.e., the cooling time) in the 
optimized case is about half that of exponential-sweep 
cooling for 60 s. The shorter sweeping time in the opti- 
mized cooling is advantageous for BEG experiments since 
the long time evaporative cooling easily becomes ineffi- 
cient due to the instability of the magnetic trap. Fur- 
thermore, it is found that the optimized frequency sweep 
becomes rapid at the final stage of the cooling process, 
i.e., after about 27 s, and such a sweep is quite opposite 
to the exponential one which slows down over time. As 
shown in the inset of Fig. 1, the optimization procedure 
accelerates the evaporation around the BEG transition 
point. Similar acceleration of evaporation at the final 
stage of optimized cooling has been also demonstrated 
by the previous studies based on kinetic theory of a clas- 
sical gas |l|,|l|. 
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FIG. 3. Cooling trajectory through phase-space for op- 
timized cooling (solid) and for exponential-sweep cooling 
(dashed). The phase-space density p is normalized by the 
critical value pc = 2.612. The BEC transition occurs when 
each trajectory passes the horizontal line where p/pc = 1- The 
arrow indicates the direction of time. Several corresponding 
time values are added to the trajectories to demonstrate the 
time evolution of the system. 



FIG. 2. Time evolution of (a) the total number of trapped 
atoms A*' and (b) temperature T, for optimized cooling (solid) 
and for exponential-sweep cooling (dashed). The dots indi- 
cate the points at which the BEC transition occurs. 

In Figs. 2(a) and 2(b), we show the time evolution of 
the number of trapped atoms N and that of temperature 
T, respectively. In the optimized cooling, both N and 
T rapidly decrease at the final stage of cooling due to 
the rapid evaporation mentioned above. On the other 
hand, TV and T decrease monotonically in the exponential 
cooling case. At the BEC transition point indicated by 
the dots in the figure, both the number of trapped atoms 
and temperature are calculated to be iV = 7.2 x 10^ and 
T — 1.2 //K at 29.01 s for the optimized cooling, and N = 
7.1 X 10^ and T = 0.56 at 48.9 s for the exponential 
cooling. 

Next, we show how the optimization affects the cool- 
ing efficiency in the whole evaporative cooling process. 
In Fig. 3, the cooling trajectories of both cases are plot- 
ted in the N-p plane and several corresponding time val- 
ues are also added to the cooling trajectories to demon- 
strate the time evolution. At the early stage of cool- 
ing where the loss due to background gas collisions is 
dominant, both trajectories are quite similar, indicating 
that the exponential sweep of rf-field frequency is suf- 
ficiently good. However, the difference becomes larger 
as the trajectories approach the quantum degenerate re- 
gion where the phase-space density exceeds the critical 
value of Pc ~ 2.612 and Bose statistics of atoms becomes 
dominant. Clearly, more atoms are lost in conventional 



exponential-sweep cooling than in the optimized case. 
For example, we have ten times larger number of atoms in 
the optimized cooling at the BEC transition point when 
the trajectory reaches the horizontal line of pj Pc — 1- 
Note that the point where the two trajectories split up 
in Fig. 3 (i.e., p/pc — 0.01) corresponds to the time of 
around 28 s where the rf-field frequency curves cross in 
Fig. 1. The sudden increase in phase-space density just 
after the BEC transition in both trajectories reflects the 
rapid growth of condensate due to bosonic stimulation. 

The number of condensed atoms, iVo, is largely in- 
creased by the optimization. In the optimized cooling, Nq 
increased monotonically and we obtained Nq = 2.1 x 10^ 
after 29.04-s cooling. In the exponential-sweep cool- 
ing, on the other hand, iVo had the maximum value of 
1.74 X 10^ at 57.0 s and then gradually decreased. We 
finally obtained A^o — 1-65 x 10^ for 60-s exponential 
sweep cooling. More than a 100-fold enhancement of the 
number of condensed atoms was achieved in our present 
calculations. The lifetime of such a large condensate is 
evaluated to be less than 100 ms due to the enormous 
three-body recombination loss, as we will see in Fig. 4. 
However, it is possible to resolve this problem by adia- 
batically expanding the trapping potential. We also add 
that the value of A'o — 1.7 x 10^ quantitatively agrees 
with the experimental results reported in Ref. ||l9[| and 
confirms the validity of our numerical calculations. 

The undesirable loss of trapped atoms during the cool- 
ing process gives us the physical origin of the rapid fre- 
quency sweep in the optimized cooling. Figure 4 shows 
the time evolution of loss rates of total number of trapped 
atoms at the final stage of the optimized cooling. It is 
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FIG. 4. Time evolution of loss rates at the final stage of 
optimized evaporative cooling, loss due to the background 
gas collision (dot-dashed), the dipolar relaxation process 
(dashed), and the three-body recombination process (solid). 



seen that the loss due to the three-body recombination 
greatly increases and selectively removes atoms with high 
density (i.e., atoms with lower energy). This means that 
the serious heating of the system is caused by three-body 
recombination loss. In order to get a large number of con- 
densed atoms, evaporative cooling should therefore be 
accelerated to compensate for such undesirable heating. 

Here, we briefly discuss the validity of the quick rether- 
malization assumption in the rapid evaporation necessary 
for optimized cooling. The evolution of the system in the 
evaporative cooling process is well scaled by the charac- 
teristic time T given by the inverse of normalized loss 
rate as r = {—Nioss/ N)~^ . The thermalization of the 
system, on the other hand, is characterized by the elastic 
collision time Tcou = {nav)~^; n is the mean atom den- 
sity, a = Stto^ is the elastic collision cross section, and 
V = AyJksT /Tim is the mean velocity of relative motion 
of the atoms . At the BEC transition point in the op- 
timized cooling, these characteristic times are calculated 
to be r = 140 ms and TcoU = 0.4 ms, respectively. The 
relation r ^ Tcoii confirms that quick rethermalization 
is still valid during the rapid evaporation that occurs for 
optimized cooling. 



efficient evaporative cooling based on general sweeping 
functions of rf-field frequency, which is useful for many 
experimentalists who want to achieve large condensates. 
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FIG. 5. (a) Sweeps of rf-field frequency consisting of two 
segments: an identical exponential sweep for 27 s and a se- 
quential linear sweep for three different times which are 1,5, 
and 10 s respectively. The exponential sweep is equal to the 
corresponding part of the dashed line in Fig. 1. The inset 
shows the linear-sweep region with an expanded scale, (b) 
Cooling trajectories through phase-space for the evaporative 
cooling using the sweeping functions in (a). The results of 
both optimized cooling and exponential cooling in Fig. 3 are 
also plotted for comparison. The phase-space density p is nor- 
malized by the critical value pc — 2.612. The BEC transition 
occurs when each trajectory passes the horizontal line where 
p/ pc ~ 1. The arrow indicates the direction of time. 



B. Efficient evaporative cooling based on general 
sweeping functions of rf-frequency 

The optimized frequency sweep in Fig. 1 is the com- 
plicated function of the experimental parameters such 
as the temperature, the density of gas, the atomic coUi- 
sional parameters, the magnetic trapping potential, and 
so on. This means that we need to calculate the opti- 
mized cooling for each BEC experiment depending on its 
experimental conditions. Here we try to demonstrate the 



From the results in Figs. 3 and 4, the exponential sweep 
is highly efficient until the three-body recombination loss 
becomes dominant. The rapid frequency sweep at the fi- 
nal stage of the optimized cooling in Fig. 1 might be 
approximated by a linear function. Thus it is expected 
that a combination of an exponential sweep of rf fre- 
quency and a rapid linear one realizes the high cooling 
efficiency. To test this cooling method, we carried out 
the time-evolution calculation by applying the same nu- 
merical parameters used in the optimization calculation 
of the proceeding subsection. As shown in Fig. 5(a), the 
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frequency is first swept exponentially from 35 MHz to 
2.26 MHz for 27 s and then swept linearly to 0.67 MHz 
for three different times of 1 s, 5 s, and 10 s respectively. 
The exponential curve in Fig. 5(a) is in conformity with 
the corresponding part of the dashed line in Fig. 1. 

Figure 5(b) shows the cooling trajectories in the iV- 
p plane for the evaporative cooling using the sweeping 
functions of Fig. 5(a), together with those for the opti- 
mized and exponential cooling plotted in Fig. 3. We can 
see that the linear frequency sweep at the final stage of 
the evaporative cooling strongly improves the cooling ef- 
ficiency as expected. Especially, the cooling trajectory 
in the case of 27-s exponential and 1-s linear frequency 
sweep (dashed line) is very close to the optimized one 
(solid line) . The final number of condensed atoms in this 
cooling method is calculated to be A'o = 1.8 x 10^ and 
comparable to the value in the optimized cooling such as 
A^o = 2.1 X 10^. One can achieve the evaporative cooling 
with sufRciently high efficiency by an initial exponential 
and a sequential rapid linear sweep of rf-field frequency. 

C. Experimental realization of efficient evaporative 
cooling 

Recently, a group at Gakushuin University performed 
the efficient evaporative cooling for ®^Rb atoms with 
F=2, mp ~ 2 trapped state j2^]. They applied a fre- 
quency sweep similar to the optimized one which we cal- 
culated according to their experimental situations. The 
time of evaporative cooling was divided into three seg- 
ments and in each segment the optimized sweeping func- 
tion was approximated by the linear one: the frequency 
was swept linearly like 32 ^13 2 ^1.1 MHz for 10, 
10, and 5 s respectively. A large number of condensed 
atoms, iVo = 1.4 x 10^, were finally obtained after this 
25-s evaporative cooling starting from the experimental 
conditions such as initial temperature T = 380 p.K, ini- 
tial total number of atoms N — 4x10^, and initial phase- 
space density p = 5.8 x 10~*. 

On the other hand, we theoretically investigated this 
evaporative cooling. It was found that the cooling tra- 
jectory through phase space is close to that of optimized 
cooling and the final 5-s linear sweep of rf-frequency en- 
hances the phase space density more than four orders 
of magnitude. The number of condensed atoms for this 
25-s evaporative cooling was calculated to be 1.7 x 10^. 
Our calculated results agree well with the experiment at 
Gakushuin University. 

IV. CONCLUSION 

We have theoretically studied the optimization of evap- 
orative cooling on the basis of quantum kinetic theory of 
a Bose gas, focusing on BEG experiments with rubidium 



atoms. The calculated results indicate that the acceler- 
ation of evaporative cooling around the BEG transition 
point is very effective against serious three-body recombi- 
nation loss. The number of condensed atoms is expected 
to be largely increased by this optimization procedure. 
Furthermore, it is shown that highly efhcient evapora- 
tive cooling can be achieved by an initial exponential 
and then a rapid linear sweep of frequency. These results 
are consistent with the experiments and provide a use- 
ful guideline to experimentalists aiming to achieve large 
condensates in ^^Rb atoms. 

In the present work, we optimized the sweeping func- 
tion of the rf-field frequency only. A more complete opti- 
mization together with the adiabatic change of trapping 
potential curvature might produce a large and long-lived 
condensate. On the other hand, our theory neglects the 
influence of gravity which reduces the cooling efficiency at 
low temperatures An interesting future work would 
be to include the gravity in our present kinetic analysis 
and study how the optimized cooling trajectory is modi- 
fied by it. 
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